Valley-Polarized Quantum Anomalous-Hall Effects in Silicene 
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We demonstrate a new quantum state of matter — valley-polarized quantum anomalous Hall effect 
in monolayer silicene. In the presence of extrinsic Rashba spin-orbit coupling and exchange field, 
monolayer silicene can host a quantum anomalous-Hall state with Chern number C = 2. We find 
that through tuning Rashba spin-orbit coupling, a topological phase transition gives rise to a valley 
polarized quantum anomalous Hall state, i.e. a quantum state which exhibits electronic properties 
of both quantum valley-Hall effect (valley Chern number Cv = 3) and quantum anomalous Hall 
effect with C = — 1. 

PACS numbers: 73.22.-f, 73.43.-f, 71.70.Ej, 85.75.-d 

Silicene, the counterpart of graphene [I] for silicon, has 
a honeycomb geometry and low-buckled structure. Due 
to its special Dirac-cone structure in low-energy spec- 
trum and promising applications in nanoelectronics [2H7J , 
silicene has attracted much attention both theoretically 
and experimentally. Distinct from graphene, silicene pos- 
sesses a stronger intrinsic spin-orbit coupling and a con- 
siderable bulk band gap can open at the Dirac points. 
Therefore, silicene becomes a good candidate to realize 
the quantum spin Hall state [7J |S] , a quantized response 
of a transverse spin current to an electric field [SHU]- 
On the other hand, a tunable extrinsic Rashba spin-orbit 
coupling from the mirror symmetry breaking about the 
silicene plane destroys this effect [HI [TDJ. Interestingly, 
this helps establish another striking topological phase — 
quantum anomalous Hall effect (QAHE) [T51419J . Unlike 
the quantum Hall effects from Landau-level quantization, 
QAHE originates from the joint effects of spin-orbit cou- 
pling and local magnetization. Although QAHE has been 
proposed for over twenty years, there was only one ev- 
idence of realization in a magnetic topological insulator 
until recently }19j . 

Similar to real spin, valleys K and K' in honeycomb 
structures provide another tunable binary degrees of 
freedom to design valleytronics. By breaking the in- 
version symmetry, e.g. the introduction of staggered 
AB sublattice potentials in honeycomb lattices, a bulk 
band gap opens to host a quantum valley-Hall effect 
(QVHE) [20 24 characterized by a valley Chern num- 
ber Cy = Ck — Ck' ■ Since both intrinsic and extrinsic 
spin-orbit couplings in silicene are comparable, one can 
expect that these effects should result in new topologi- 
cal phenomena. It has turned out that manipulating the 
spin-orbit couplings is one of the most effective means to 
control the properties of silicene. For example, an inter- 
esting valley-polarized metallic state has been reported 



in silicene recently [25] . Because the intrinsic and Rashba 
spin-orbit couplings give different responses on the con- 
duction and valence bands of silicene, it is natural to 
imagine that a topological phase including both proper- 
ties of QAHE and QVHE may exist in silicene. 

In this Letter, we report a valley-polarized quantum 
anomalous Hall phase in monolayer silicene. When the 
time reversal symmetry is broken from the exchange field, 
we find that the combination of intrinsic and Rashba 
spin-orbit couplings can result in a new topological phase, 
i.e., a bulk band gap closing and reopening occurs along 
with the increasing of Rashba spin-orbit coupling at fixed 
intrinsic spin-orbit coupling and exchange field. Through 
analyzing the resulting Berry curvatures for the occupied 
valence bands, we show that the nonzero Chern number 
indicates a quantum anomalous-Hall phase. Surprisingly, 
we further find that valleys K and K' give rise to differ- 
ent Chern numbers, i.e. Ck = 1 and Ck> = — 2. This 
imbalance signals a quantum valley-Hall phase with val- 
ley Chern number Cv = 3. Our finding not only provides 
a platform to design low-power electronics but also ad- 
vances the application of silicene-based valleytronics. 

In the tight-binding approximation, the Hamiltonian 
for silicene in the presence of spin-orbit couplings and 
exchange field can be written as [7J [5] : 
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where c ia {ci a ) is a creation (annihilation) operator for an 
electron with spin a on site i. The first term represents 
the nearest neighbor hopping term with hopping energy 
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FIG. 1: (Color online) Evolution of band structures of the bulk [panel (a)-panel (f)] and zigzag-terminated [panel (g)-panel 
(1)] silicene as a function of Rashba spin-orbit coupling ta at fixed intrinsic spin-orbit coupling tso and exchange field M. (a) 
tR — 0. Bulk energy gaps open around the K and K' Dirac points. The size of bulk gap near valley K is exactly the same 
as that near valley K' . Panels (b)-(f): tR — 0.01,0.044,0.06,0.08,0.12, respectively. Along with the increasing of tR, the bulk 
gap around valley K gradually increases, while the bulk gap near valley K' closes twice [see panels (c) and (e)] and reopens 
twice [see panels (d) and (f)]. In panels (g)-(l), the valley-associated gapless edge modes at valley K is unchanged, but those 
for valley K' changes, i.e. there are two/one pairs edge modes after the bulk gap reopens. Other parameters are set to be 
to/t = 0.07, tso/t = 0.10, and M/t = 0.50. 



t. The second term is the effective spin-orbit coupling 
involving the next-nearest neighbor hopping with ampli- 
tude to- Vij — di x dj/\di x dj\ , where di and dj are two 
nearest bonds connecting the next-nearest neighbor sites. 
The summation over (...) (((...))) runs over all the near- 
est (next-nearest) neighbor sites. The third and fourth 
terms are respectively the intrinsic and Rahsba spin-orbit 
couplings with tso an d tji the corresponding spin-orbit 
coupling strengths, dij = dij/\dij\, where d\j represents 
a vector from site j to i and \Xij = ±1 for A or B site. 
The last term is an exchange field M, which arises from 
the interaction with a magnetic substrate. 

By transforming the real-space Hamiltonian in Eq. (fTl) 
into a 4 x 4 matrix H(k) for each crystal momentum k on 
the basis of {ipAf, Vvu, ipRf> ipBi}, the band structure of 
bulk silicene can be numerically obtained by diagonaliz- 
ing H(k). Panel (a)-panel (f) in Fig.lllplot the evolution 
of the bulk band structure along with the increasing of 
Rashba spin-orbit coupling tR. Other parameters are set 
to be t /t = 0.07, tso/t = 0.10 and M/t = 0.50. Thanks 
to the absence of intervalley scattering, valleys K and 
K 1 remain good quantum numbers and distinguishable. 
In the absence of Rashba spin-orbit coupling tn/t = 0, 
one can see that bulk band gaps open at valleys K and 
K' [55]. When the Fermi level lies inside the bulk en- 
ergy gap, it is known that this insulating state is a quan- 
tum anomalous Hall insulator, which is characterized by 
a nonzero Chern number C 26-28 calculated from 
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where £l n is the momentum-space Berry curvature for the 



nth band [13 \M I 
Q n (k) =-J2 
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The summation is over all occupied valence bands in the 
first brilloin zone blow the bulk energy gap, and v x r y ) is 
the velocity operator along x(y) direction. The absolute 
value of C corresponds to the number of gapless chiral 
edge states along any edge of the two-dimensional sys- 
tem. When tn/t = 0, the Chern number obtained from 
the tight-binding Hamiltonian is C — 2. Furthermore, by 
calculating the Chern number contribution of each val- 
ley based on a continuum model Hamiltonian, it is found 
that Ck = Ck> =1- In this case, the corresponding valley 
Chern number Cy = Ck — Ck' is vanishing. As plotted 
in panels (b)-(f), when the Rashba spin-orbit coupling 
gradually increases from t^/t = 0.01 to tn/t = 0.12, we 
find that the bulk band gap at valley K always increases. 
To our surprise, the corresponding bulk band gap at val- 
ley K' closes and reopens twice, signaling two possible 
topological phase transitions. 

Let us now explore the topological properties of the 
resulting two phases based on the Chern number calcu- 
lation described in Eq. pi). In the first phase after the 
gap reopening shown in Fig. md), the Chern number is 
obtained to be C = — 1. Since the bulk gap near valley 
K does not close, the resulting Chern number keeps un- 
changed. However, the corresponding contribution from 
valley K' becomes Ck 1 = — 2. Naturally, this imbalance 
of Chern number contribution between valleys K and K' 
gives rise to a nonzero valley Chern number Cy = 3. This 
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FIG. 2: (Color online) (a) Contour of Berry curvature distribution in (k x , k y ) plane for the valley-polarized QAHE. (b) Berry 
curvature distribution as a function of k x at fixed k y = 2-7r/\/3. (c) Valley-associated edge modes for the valley-polarized 
QAHE. 



means that the new topological insulating phase is both a 
QAHE and a QVHE. Such a state is robust against any- 
kind of impurities (magnetic or not) and carries a net 
valley-polarized current. For simplicity, hereinafter we 
name it as a valley-polarized QAHE. To our knowledge, 
this is the first time such phase is reported. To further 
explore this new phase, in Fig. u[a) we show the contour 
of its Berry curvature distribution in the (/c x , k y ) plane. 
And Fig. [2] (b) plots the Berry curvature as a function 
of k x at fixed k y = 27r/v3. One can find that the Berry 
curvatures are mainly localized around the valley points, 
and obviously the Berry curvature density near valley K 
is different from that near valley K'. This is a direct con- 
sequence of the inequality of the Chern number carried 
by valley K and K'. 

In another new phase after the second topological 
phase transition for even larger ta [see in Fig. [TTf )] , the 
Chern number becomes C — 2 again, where Ck' = 1 for 
valley K', indicating that this phase is the same as the 
phase without Rashba spin-orbit coupling. In this way, 
we have realized a Chern-number-tunable QAHE by con- 
trolling the Rashba spin-orbit coupling strength. In ad- 
dition to the Chern number, gapless edge mode inside 
the bulk energy gap provides a more intuitional picture 
to characterize the QAHE. 

Figures fljg)-[l|l) plot the one-dimensional band struc- 
ture of zigzag-terminated silicene ribbon. One can ob- 
serve that near valley K there is always one pair of edge 
mode for any Rashba spin-orbit coupling strength. While 
at valley K' the corresponding edge modes vary along 
with the topological phase transitions, i.e., after the first 
and second topological phase transitions there are respec- 
tively two and one pair of edge modes. Through studying 
the wavefunction distribution of the gapless edge states 
inside the bulk band gap shown in Fig. flTj) and from 
the energy dispersion, we find that the edge modes of 
the valley-polarized QAHE have the form as plotted in 
Fig. [2] (c): (1) there are three edge states localized at 



each boundary; (2) For the upper (lower) boundary, two 
edge states encoded with valley K'(K) propagate from 
left (right) to right (left) while one encoded with valley 
K (K ! ) counterpropagate from right (left) to left (right); 
(3) a net valley-polarized current associate with valley 
K' (K) flows to the right (left) direction. 

In below, we present a simple physical understanding 
of the formation of this valley-polarized QAHE using the 
spin texture picture. By analyzing the spin textures of 
the real-spin s and the pseudo-spin textures of AB sublat- 
tices <r [3T] , we show that the topological charges carried 
by the real-spin or AB sublattice pseudo-spin for valleys 
K and K' are 
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where '1' and '2' labels the two valence bands below the 
bulk band gap. From the above finding, one can conclude 
that the net Chern number is mainly contributed from 
the real-spin texture carried by the 2 nd valence band. 

Finally, we provide a phase diagram in the (tso, til) 
plane. Figure [3] clearly shows that there are three topo- 
logical phases separated by two gray lines representing 
the topological phase transition. Here, we use color to 
signify the size of the corresponding bulk band gap. From 
the Chern number calculation, we show that the Chern 
number in phase-I is identical to that in phase-Ill, i.e., 
C = 2. However, in phase-II the Chern number is C = — I 
and the valley Chern number is Cy = 3, which is exactly 
the newly found valley-polarized QAHE. 

In summary, in the monolayer silicene we have found 
by tuning the Rashba spin-orbit coupling, a new topologi- 
cal phase — the valley-polarized quantum anomalous Hall 
state. Different from a conventional quantum anoma- 
lous Hall state, the valley-polarized quantum anomalous 
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FIG. 3: (Color online) Phase diagram in the (tso, tii) plane. 
Other parameters are set to be to/t — 0.07 and M/t = 0.50. 
The grey lines separate three topological phases for the single 
layer silicene, and colors are used to indicate the size of the 
bulk band gap. In phase-I and phase-Ill, the Chern numbers 
are identical C — 2. While in phase-II, the Chern number is 
C = -1 and Cy = 3. 



Hall state has a quantized Chern number C = — 1 and 
a nonzero valley Chern number Cy = 3. Therefore, the 
new topological phase possesses the properties of both 
quantum anomalous Hall effect and quantum valley Hall 
effect, which is a good candidate to design dissipationless 
valley tronics. The Rashba spin-orbit coupling can serve 
as a topological switch to drive the monolayer silicene 
from a conventional quantum anomalous Hall phase to 
a valley-polarized quantum anomalous Hall phase, which 
can be realized by controlling the adatom coverage in 
monolayer silicene. 
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